function OBSplot(logFlag)

%% observation data
obs{1}=[8 12.6 28.6 9.8 22.4 2.1 2.3 5.1 8.5];%AM0644x3	
obs{2}=[8.8 7.1 65.1];%Arp148x1	
obs{3}=[2.7 9.8 2.9 1.9 3.8 1.6 4.3 1.4 1.2 2.8 7.5 1.4 2.2 1.6];%Arp143
obs{4}=[3.36 11.86 11.15 65.49 7.43 6.37 17.7 2.66 15.93 19.47 2.3 1.19 2.47 1.45 0.86];
%Cartwheel
obs{5}=[0.69 3.92 24.49 1.2 1.03 0.53 3.14 0.99 3.96 8.77 0.89 0.65];
%N922
obs{6}=[7.52 3.39 1.47 2.06 3.61 2.43 5.63 6.25 5.54];%Arp 147-Ring
obs{7}=[64.1 0.5 0.6 0.7 0.8 0.6 0.3 0.8 24.1];%Arp 284
obsall=[];
SFRgalaxy= [2.6 2.5 2.3 20 8.0 4.1 4];

%% analysis and plot
% for i=1:7
%     obsall=[obsall; [obs{i} ; ones(size(obs{i}))/SFRgalaxy(i)  ]' ];
% end
% obsall(:,1) =log10( obsall(:,1) ) + 39;
% obssort=sortrows(obsall,1);
% % num=linspace(length(obssort),1,length(obssort));
% num= cumsum( obssort(:,2)','reverse' ) ./7;
% % num=log10(num);
% 
% % num = num * 43.5;
% if(logFlag)
%     num=log10(num);
% end
% plot(obssort(:,1),num,'ro')



for i=1:7
    obsall=[obsall,obs{i} ];
end
obsall =log10( obsall) + 39;
obssort=sort(obsall);
num=linspace(length(obssort),1,length(obssort));
% num=log10(num);
if(logFlag)
    num=log10(num);
end
plot(obssort,num,'ro')



%% setplot
% opt.LineStyle = {''}; % three line styles
% opt.Markers = { 's'};
% % opt.MarkerSpacing = [15, 15, 15];
% opt.Legend = {'Observation', '\theta = 45^o', '\theta = 90^o'}; % legends
% 
% % Save? comment the following line if you do not want to save
% opt.FileName = 'plotMarkers.png'; 
% 
% % create the plot
% setPlotProp(opt);

end